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Abstract 

We use an improved nonperturbative variational method to investigate the 
phase transition of the Gross-Neveu model. It is shown that the variational 
procedure can be generalized to the finite temperature case. The large N 
result for the phase transition is correctly reproduced. 
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I. INTRODUCTION 


A survey of the literature concerning analytical nonperturbative methods shows a grow¬ 
ing interest in variational methods such as the gaussian approximation and the optimized 
linear 6 expansion [Q. Recently, an alternative variational method has been improved 
and successfully tested in the Gross-Neveu (GN) model both in the large and hnite 
N limits. Basically this new approximation differs from the existing optimized linear 6 ex¬ 
pansion in two aspects: it provides a clear way of carrying on the renormalization program of 
the theory and improves the calculation with a simple but powerful analytical method which 
allows to attain arbitrarily large and even inhnite order of perturbation in the parameter 
S. Our goal is to show that this calculational scheme, especially the variational procedure, 
still works when external parameters such as the temperature are considered. From the 
renormalization point of view the temperature does not create any new difficulties and the 
elimination of divergencies can be carried out exactly as at zero temperature ^j. However, 
previous experience with the related S expansion shows that the application of variational 
methods to the study of phase transitions is not necessarily straightforward 0. 

Here, our aim is to check the ability of the variational procedure used in Ref. in 
generating accurate results for the phase transition, especially near the critical temperature. 
In this hrst hnite temperature application we restrict ourselves to the large N limit of the 
GN model where the calculational scheme can be set up more clearly. For comparison 
we consider the results of Refs. §] and as being “exact” in the large N limit, despite 
their disagreement with existing theorems [0 for symmetry breaking/restoration in 1-1-1 
dimensions. This interesting problem, which has been studied by several authors using 
different models [11-15], will be explicitly treated in an extension of the present work fT^ . 

In the next section we review the usual large N result for the temperature dependent 
fermionic mass and perform the variational calculation to lowest orders. We shall see that, 
as in the zero temperature case, the variational procedure fails. In Section HI we follow Ref. 
[^ to rectify the situation by performing the variational calculation to all orders. When this 
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is done the large N result is completely reproduced. The conclusions and future perspectives 
are presented in section IV. 


II. THE VARIATIONAL CALCULATION TO LOWEST ORDERS 


The variational calculation starts with the addition of an arbitrary bare mass (mo) to 
the original massless Gross-Neveu 


N 


N 


i=l i=l ^ \i=l / 


( 1 ) 


where go is the bare coupling constant (in the following we shall suppress the summation 
over the index i). The relation to the linear 6 expansion and other variational methods 
becomes clear by performing the substitutions 


mo ^ mo(l - S) 


( 2 ) 


2 r 2 

9o ^9o • 


(3) 


These will be done at a later stage in order to avoid the explicit evaluation of Feynman 
graphs which differ only by 6mo insertions. To perform hnite temperature calculations in 
the imaginary time formalism one does the following substitutions 


dpo 

2tt 




(4) 


PO ^ ? 


(5) 


where for fermions 


uJn = T{2n + l)7r . 


( 6 ) 


The sum over Matsubara’s frequencies can be performed with 


T ^ ln(a;^ + ml) = E + 2Tln[l + exp{-E/T)] + c , 


(7) 
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where c is a ii^-independent constant. In 1+1 dimensions 


E = (pi+mo)5 . 


The remaining space integral is evaluated in 1 — e dimensions using conventional dimensional 
regularization techniques. 

In the large N limit a perturbative calculation of the fermionic mass {Mp) to OIqq) 
yields 


= mP + glmP + g^mP + O(^o) > 


where mP = mo, 


( 1 ) N 


my = 


ml " [(47r) 2 r(e/2) - 4/f (yo)] , 


mP = —ml-^p{47ryr\e/2){l-e) 

+4(47r)tr(e/2) [(3e - 2)Ipyo) + (mo/T)2(l - e){Ipyo) + Ipyo)) 
+16/f (yo) [if (yo) - {mo/Tf{I^(yo) + ij (yo))] } • 


The temperature dependent integrals are 


ipyo) = 


roo 

-{27iyoy / dx 
Jo 


r(l/ 2 -e/ 2 )^ Jo ^ y2y2[i y exp{x‘^ + y^)2] 


il (yo) = 


roo 

-{2nyoY / dx 
Jo 


r(l/ 2 -e/ 2 ) Jo ^2^2 [1 + exp(a ;2 + 1/^)2] 


4 (Vo) = 


r(l/2-e/2) 


roo 

{2nyoy / dx 
Jo 


X exp(a;^ + I/q) 2 

+ 1 / 0 ) [1 + exp(a;2 + ylpf 


with yo = mo/T and x = pi/T. These integrals are related to each other via 


Myy^iT) 

dyi 


_ o 6 ^140 + ) _ jT ,jT 

€, 9. ~ -'2 -'s 
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The nonperturbative large N calculation consists in the evaluation of all cactus diagrams 


which can be summed up as 


Mf = mo [l- ^glMr [( 47 r)tr(e/ 2 ) - 4Jf (|/^)] | , 

where yp = Mp/T. The bare and renormalized parameters are related via 


(16) 


mo = Z„,m , 


( 17 ) 


and 


9o = , 


(18) 


where /i is the arbitrary scale introduced by dimensional regularization. The renormalization 
constants are 

-1 


Zm — Zg — 


1 + 


g^N 


vre 


(19) 


Substituting Eqs. ([T7|) and (^) into Eq. ([I6D yields the hnite expression for the dimensionless 
quantity Mp/Aj^ 

-1 




m 


N 


^MS ^MS 1 ^ 






+ 2Ii(yp) 


( 20 ) 


where Ajj^ = y,expl—ir/{Ng'^)] and y = p.( 47 r )2 exp(— 7 ^/ 2 ) with 7 ^ = 0.577215. Equation 
(^) satishes the RG equation 

d ( d d d \ 

^ Mp=iy— + gl3{g)—-m'ym{g)^]Mp = 0 , ( 21 ) 

ay \ ay og am J 

where 

Pig) = -^g"^ and jrnig) = —g"^ ■ (22) 

/TT 71 

Then, in the chiral limit (m = 0), the hnite temperature mass gap equation is given by 


M^(T)=M^(0)exp[-2Jn2/F)] , 


(23) 
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where Mp{Q) = is the renormalized fermion mass, generated through dynamical sym¬ 
metry breaking at T = 0. This is just the large N result which reproduces the phase diagram 
for the fermion mass [Q. 

The variational calculation starts with the substitutions Eqs. (^) and . With this 
procedure Eq. interpolates between the original massless GN model (at 5 = 1) and a 
massive free theory (at 5 = 0). The perturbative calculation is now done in powers of the 
bookkeeping parameter 6 and extremized with respect to mo at 5 = 1 (Principle of Minimal 
Sensitivity J^). Starting with 0(5) one has at 5 = 1 


M^p\mo) = mp’ + gorrip’ - mo 




2_(i) 


dm^P 2 ( 1 ) 


dmn 


= %rny 


(24) 


which has no nontrivial extremum in mo at T = 0 nor at hnite temperatures. At 5 = 1 the 
0(5^) fermionic mass is given by 




(25) 


The extremization with respect to mo in the limit e —0 does not give any useful information 
since the term g^m^P, which has no nontrivial extremum, dominates. This behaviour, which 
persists to higher orders, has also been noted at zero temperature and we hnd that the 
situation does not change at hnite temperatures. 


III. THE VARIATIONAL CALCULATION TO ALL ORDERS. 

In this section we shall follow Ref. to perform an all order variational calculation in the 
large N limit eliminating the optimization problem encountered in the previous section. The 
general philosophy within variational methods is to start with a trial value which is expected 
to be reasonably close to the true value of the physical parameters. In our case this means 
that we can start by formulating a nonperturbative ansatz which already resums a good 
part of the RG behaviour of the fermionic mass before launching into the actual variational 
calculation. Of course this is a rather easy task within the large N limit, where the exact 
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answer Eq. (^) constitutes the natural choice. Performing the substitutions Eqs. (0) and 
(|) in Eq. ([T6|) we get 


MfU) = 


mo(l - 6 ) 

m 


(26) 


where we have dehned 


N 


nS) = 1 - - S)-’r [(4^)>r(£/2) - A!‘ (j/H 


(27) 


with y'p = MF{f)/T. It is now possible to perform an expansion in powers of 6 to order-n 
around the free theory (5 = 0). Using contour integration one obtains Mp to order of 
perturbation theory 


T^An), \ 1 / 7 /I ^ N"io(l —2;) 


f{z) 


(28) 


which, at 5 = 1, gives 




dz rriQ 


(29) 


^n+l 

Once Zm and Zg are applied to the bare mo and one gets the hnite expression for the 
dimensionless quantity Mf/A-jj^ 


MP{m) 

~A 


MS 


dz m 


2Fi I A' 


MS 


N 2 

1 + -g^z 

TT 




n ^ -1 


(30) 


As noted in the case of the anharmonic oscillator , it is possible to extract more structure 
from the limit of inhnite order by rescalling m with the order n. After distortion of the 
contour it is clear that only the vicinity oi z = 1 survives in the limit n ^ oo, which can be 
analyzed by changing variables 


l-z = 


n 


MF(m") 1 / dve^m” 

Mrs 277! t K{v) 


(31) 


Rescalling m by introducing m' = m/n we get, in the n —>■ oo limit 
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( 32 ) 


















where the integration runs counterclockwise around the negative real axis. The function 
K{v) is given by 


K{v) = In 


m"v 


K{v) 


+ 21 ^ 


m"v 


tK{v) 


(33) 


where 


m" = 


m 


'Ng^ 


-1 


A 


MS 


TT 


(34) 


with t = T/Aj^. Equations (|3^) and (|3^) , which summarize our variational approach, 
should be understood as follows: for a given variational parameter m", and a given tem¬ 
perature T, Eq. enables one to determine K{v) self consistently. The variational result 
is then given in an explicitly RG invariant way by extremizing Eq. (|3^ with respect to 
m". In general this program has to be achieved numerically. However, before doing that 
one can use the fact that the large N limit is free from infra red divergences to perform an 
analytical exploitation of the m" —> 0 limit where the integral is dominated by the v ~ 0 
region. Simple considerations show that for T < 



m"v 


^MS 

Mf{T) 


(35) 


where Mf{T) is given in Eq. Moreover, K{v) has a cut starting at a negative value of v 
and lying along the negative real axis. Hence, the integral Eq. converges exponentially to 
the expected result as m" ^ 0. As T —>■ T^, the branching point approaches the value n = 0 
merging to it at T = T^. At this point the integral becomes divergent at m" = 0 and does 
not allow more extrema for T > T^. Numerical results obtained at different temperatures 
(see Fig. 1) indicate that m" = 0 is in fact the only real extremum which leads to a smooth 
(second order) phase transition occuring at the critical temperature 


G = 0.57 . 


( 36 ) 


The standard large N structure of the phase transition is then exactly reproduced. 
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IV. CONCLUSIONS 


In this letter, we have shown that a recently proposed variational method can be snc- 
cessfnlly generalized to the hnite temperatnre domain. For comparison pnrposes, and to 
introduce the method in the study of phase transitions we have chosen to start with the 
large N limit of the Gross-Neveu model. We have seen that optimization problems encoun¬ 
tered in the zero temperature low order variational calculation persist at hnite temperatures. 
Then, by applying the all order variational calculation scheme developed in Ref. [Q we were 
able to recover the usual large N result for the phase transitionof the GN model, showing 
that the whole calculational scheme works remarkably well at all temperatures. The hnite 
N case, which is interesting due to some important theorems related to phase transitions 
in one space dimension is also technically much more complex and will be discussed in a 
forthcoming work. 
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FIGURES 

FIG. 1. Mp{T)/Mp{Q) as a function of the arbitrary parameter m" for different temperatures. 
From top to bottom the curves represent t = 0.1, t = 0.5 and t = 0.55 respectively. 
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